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Abstract 

00 . 

, We have employed the framework of Bethe-Salpeter equation under covariant instantaneous 

ansatz to calculate leptonic decay constants of unequal mass pseudoscalar mesons like Ti^, K, D, D s 
and B and radiative decay constants of neutral pseudoscalar mesons like 7r° and r] c in two photons. 
In the Dirac structure of hadronic Bethe-Salpeter wave function, the covariants are incorporated 
from their complete set in accordance with a recently proposed power counting rule. The decay 
constants are calculated with the incorporation of both Leading order and Next-to-leading order 
Dirac covariants. The results validate the power counting rule which provides a practical means of 
incorporating Dirac covariants in the Bethe-Salpeter wave function for a hadron. 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD) is the theory to describe strong interactions. How- 
ever, the large gauge coupling at low energies (long distances) destroys the perturbative ex- 
pansion. As a result, many non-perturbative approaches have been proposed to deal with this 
long distance properties of QCD, such as QCD sum rules, Lattice QCD, dynamical-equation- 
based approaches like Schwinger-Dyson equation and Bethe-Salpeter equation (BSE), and 
potential models. Since the task of calculating hadron structures from QCD itself is very 
difficult, as can be seen from various Lattice QCD approaches, one generally relies on specific 
models to gain some understanding of QCD at low energies. BSE is a conventional approach 
in dealing with relativistic bound state problems. From the solutions we can obtain useful 
information about the inner structure of hadrons, which is also crucial in treating high en- 
ergy hadronic scatterings. The BSE framework is firmly rooted in field theory, and provides 
a realistic description for analyzing hadrons as composite objects. Despite its drawback of 
having to input model-dependent kernel, these studies have become an interesting topic in 
recent years, since calculations have shown that BSE framework using phenomenological 
potentials can give satisfactory results on more and more data being accumulated. 

In this paper we study leptonic decays of pseudoscalar mesons (P-mesons) such as n, 
K, D, Ds and B, which proceed through the coupling of quark-antiquark loop to the 
axial vector current and also the two-photon decays of neutral pseudoscalar mesons such 
as 7r°and r] c which proceed through the famous quark-triangle diagrams. We employ QCD 
motivated BSE under Covariant Instantaneous Ansatz (CIA) in this paper [l|, CIA 
is a Lorentz-invariant generalization of Instantaneous Ansatz. For a qq system, the CIA 
formulation ensures an exact interconnection between 3D and 4D forms of BSE {2], S]. The 
3D form of BSE serves for making contact with the mass spectrum, whereas the 4D form 
provides the Hadron-quark vertex function for evaluation of various hadronic transition 
amplitudes through quark loop diagrams. In these studies one of the main ingredients is the 
Dirac structure of the Bethe-Salpeter wave function (BSW). The copious Dirac structure of 
BSW was already studiedd by Llewllyn Smith J] much earlier. Recent studies have 
revealed that various covariant structures in BSWs of various hadrons is necessary to obtain 
quantitatively accurate observables. It has been further noticed that all covariants do not 
contribute equally for calculation of meson observables. So it is interesting to investigate 
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how to arrange these covariants. In a recent work [2j, we developed a power counting rule 
for incorporating various Dirac structures in BSW, order-by-order in powers of inverse of 
meson mass. We have outlined the Dirac covaiants and expanded the coefficients to the 
leading order (LO), and calculated the leptonic decay constants of vector mesons (p, u, <fi, 
ijj) 2I as well as pseudoscalar mesons (tt, K, D, D s and B) j^] at this order. The results 
agree with data well. 

However, common to all the perturbative theories, it is better to calculate the next 
order(s) to the leading one and make sure it is (they are) really smaller w.r.t. the LO, before 
claiming the validation of the perturbation. At the same time, as more and preciser data 
accummulated, it is useful to arrange more available parameters inherent in our framework 
to accommodate better fits to gain more precise information of the structure of hadron. So 
the study of next-to-leading order (NLO) is natural and essential. For all the mesons, the 
pseudoscalar is the simplest in Dirac structure. As the first step, we collect the data of 
leptonic decay constants /p's for pseudoscalar mesons (tt, K, D, Ds and B), to fit three 
parameters B\ s in our framework at NLO. We found: a) NLO works better than LO. b) NLO 
corrections are smaller than those of LO (it is exceptional for its small mass, to be discussed 
later in this paper). Then with the fitted parameters we calculate the radiative decay 
constants Fp of neutral pseudoscalar mesons, ir° and r] c at NLO. We also found satisfying 
agreement with data, and fair improvement w.r.t. LO. Thus the fact that three parameters 
can give a good fit not only for 5 different cases of fp, but also giving satisfactory results for 
two cases of Fp, demonstrates the validity and robustness of this framework. These results 

n 

indicate that our power counting scheme [2J provides a practical means of incorporating 
various Dirac structures from their complete set into the BS wave function. 

In what follows, we give a detailed discussion of the fit and calculation at NLO, after a 
brief review of our framework. The paper is organized as follows: In section [Til we discuss 
the structure of BS wave function for P-mesons in BSE under CIA using the power counting 
rule. In section IIIIl we introduce the fitting to fp for pseudoscalar mesons. The radiative 
decay constants Fp for 7r°and r\ c mesons are calculated in section 4, while we conclude with 
Discussion in section [V] 
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II. THE BSW UNDER CIA 



A. BSE under CIA 

We first outline the BSE framework under CIA. We have employed for the case of scalar 
quarks for simplicity. For a qq system with an effective kernel K and 4D wave function 
$(P, q), the 4D BSE takes the form, 

i (2vr) 4 A 1 A 2 $(P, q) = J d A qK{q, q')$(P, q r ), (1) 

where A 12 = ml 2 + pf 2 are the inverse propagators of two scalar quarks, and m\ t 2 are 
(effective) constituent masses of quarks. The 4-momenta of the quark and anti-quark, 
are related to the internal 4-momentum q^ and total momentum P^ of hadron of mass M as 

Pl,2ft = ™>l,2Pn ± (2) 

where = [1 ± (m\ — m\)/M 2 ]/2 are the Wightman-Garding (WG) definitions of masses 
of individual quarks. Now it is convenient to express the internal momentum of the hadron 
q as the sum of two parts, the transverse component, q^ = g M — ^ which is orthogonal 
to total hadron momentum P (ie. q.P = regardless of whether the individual quarks are 
on-shell or off-shell), and the longitudinal component, oP^ = (q ■ P/P 2 )P fl , which is parallel 
to P. We now use an Ansatz on the BS kernel K in Eq. ([T]) which is assumed to depend on 
the 3D variables q^, q'^ p| i.e. 

K(q,q') = K(q,q'), (3) 

A similar form of the BS kernel was also earlier suggested in ref. js]). Hence, the longitudinal 
component, aP^ of q^, does not appear in the form K(q, q') of the kernel. For reducing Eq.([T]) 
to the 3D form, we define a 3D wave function 0(g) as 

+oo 

<P(q) = [ Mda§{P,q). (4) 



4 



Substituting Eq. (jlj) in Eq. ([T]), with definition of kernel in Eq. ([3]), we get a covariant 
version of Salpeter equation, 

(2nfD(q)<j ) (q) = f d 3 q' K(q, q')<J)(q'), (5) 
where D(q) is the 3D denominator function defined by 



+00 

1 f Mda 



D(q) 2m J AiA 2 ' 



(6) 



-co 



whose value can be easily worked out by contour integration by noting positions of poles in 



the complex a-plane (shown in detail in 



as, 



= 1 1 > ^1,2 = ^1,2 + 9 • (7) 

2a>i 2u; 2 

We can see that RHS of Eq. ([5]) is identical to RHS of Eq. ([T|) by virtue of Equations ([3]) 
and (jl]). We thus have an exact interconnection between 3D wave function <p(q) and 4D 
wave function $(P, q): 

A A *(P,,) = «.rffl. (8) 

lm 

We also get the Hqq vertex function T(q) under CIA for case of scalar quarks. Further 
in the process, an exact interconnection between 3D and 4D BSE [7] is thus brought out 
where the 3D form serves for making contact with the mass spectrum of hadrons, whereas 
the 4D form provides the vertex Hqq function T(q) which satisfies a 4D BSE with a natural 
off-shell extension over the entire 4D space (due to the positive definiteness of the quantity 
q 2 = q 2 — (q ■ P) 2 / P 2 throughout the entire 4D space) and thus provides a fully Lorentz- 
invariant basis for evaluation of various transition amplitudes through various quark loop 
diagrams. 
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B. Dirac structure of Hadron-quark vertex function for P-mesons in BSE with 
power counting scheme 



To obtain the form of Hadron-quark vertex function for the case of fermionic quarks 
constituting a particular meson, we first replace the scalar propagators A^ 1 in Eq. (j7j) by 
the proper fermionic propagators Sf- The Hqq vertex function T(q) now is a 4 x 4 matrix 
in spinor space. For incorporation of the relevant Dirac structures in T(q), we make use of 
the power counting rule we developed in 2], order- by-order in powers of inverse of meson 
mass [2]. Our aim of developing the power counting rule was to find a "criterion" so as to 
systematically choose among various Dirac covariants from their complete set to write wave 
functions for different mesons (vector mesons, pseudoscalar mesons etc.). 

As far as a pseudoscalar meson is concerned, its hadron-quark vertex function which 
has a certain dimensionality of mass can be expressed as a linear combination of four Dirac 



covariants J4], each multiplying a Lorentz scalar amplitude, as function ofq-P. We note that 
in the expression for CIA vertex function in equation (E3), the factor D(q)<f)(q) is nothing but 
the Lorentz- invariant momentum dependent scalar which depends on q 2 , P 2 and q ■ P and 
has a certain dimensionality of mass. However the Lorentz-scalar amplitudes multiplying 
various Dirac structures in [5| have different dimensionalities of mass. For adapting this 
decomposition to write the structure of Hqq vertex function T(q) for a particular meson, we 
re-express this function by making these scalar amplitudes dimensionless by weighing each 
covariant by an appropriate power M, the meson mass. Thus each term in the expansion 
of r(g) is associated with a certain power of M and hence in detail we can express the 
hadron-quark vertex, T(q) as a polynomial in various powers of 1/M: 

T p (q)=n p ^- i N P D(q)cf>(q), (9) 

with 

n p = l5 B - z 7 s(7 • P)^ - t7fe(7 • ?)§ - 7fe[(7 ■ P)(l ■ q) ~ (l ■ q){l ■ (10) 

where Bi (i = 0, ...3) are four dimensionless coefficients to be determined. Since we use 
constituent quark masses, where quark mass m is approximately half of the hadron mass 
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M, we can use the ansatz 



q«P~M (11) 

in the rest frame of the hadron (however we wish to mention that among all the pseudoscalar 
mesons, pion enjoys the special status in view of its unusually small mass (M < Aqcd) and 
its case should be considered separately). Then each of the four terms in Eq. (l9][T0l) would 
again receive suppression by different powers of 1/M. Thus we can arrange these terms 
as an expansion in powers of 0(1/M). We can then see in the expansion of Q p , that the 
structures associated with the coefficients Bq, B\ have magnitudes 0(1/M°) and are of 
leading order, while those with B 2 , B 3 are O^i/M 1 ) and are next-to- leading-order. This 
naive power counting rule suggests that the maximum contribution to the calculation of any 
pseudoscalar meson observable should come from the Dirac structures 75 and 275 (7 • P) /M 
associated with the constant coefficients B and B\ respectively, followed by the other two 
higher order covariants associated with coefficients B 2 and B 3 . In general, the coefficients Bi 
of the Dirac structures could be functions of q-P, and hence can be written as a Taylor series 
in powers of q ■ P. However the coefficients used here are dimensionless on lines of [2|. So 
they are in fact function of q ■ P/M 2 . Then the leading order contribution of the coefficients 
are the case when the B^s are constant. In this paper, we assume the coefficients are smooth 
functions of q ■ P/M 2 , so to NLO, we only consider the terms of eq.(10), with the coefficients 
Bi constant. Because the normalization of the BSW can be fixed (see below), B Q here can 
be taken to be 1. So we totally have 3 parameters to be fitted at NLO, comparing to one 
parameter at LO. In a similar manner one can express the full hadron-quark vertex function 
for a scalar and axial vector meson also in BSE under CIA. At the same time, the restriction 
by charge parity on wave function of eigenstate should also be respected. Further, to get 
the complete set of the Dirac structures for a certain kind of meson, the restriction by the 
(space) Parity have been employed; and it is easy to see that the requirements of the space 
^arity and the charge Parity are the same for the vertex as well as the full wavefunction 
lp| . In this work to calculate the leptonic and radiative decay constants, we take the form 
of hadron-quark vertex as in Eqs. (9) and (10) which incorporates LO as well as NLO 
covariants and see the relative importance of various covariants. 
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C. BSE Kernel and the scalar wave function 



From the above analysis of the structure of vertex function Hqq, we notice that the 
structure of 3D wave function <f)(q) as well as the form of the 3D BSE are left untouched and 
have the same form as in our previous works which justifies the usage of the same form of the 
input kernel we used earlier [2j . Now we briefly mention some features of the BS formulation 



employed. The structure of BSE is characterized by a sing 



e effective kernel arising out 



of a four-fermion laerangian in the Nambu-Jonalasino 11, 12 ] sense. The formalism is 

n 

fully consistent with Nambu-Jona-Lasino [11( picture of chiral symmetry breaking but is 
additionally Lorentz- invariant because of the unique properties of the quantity q 2 , which is 
positive definite throughout the entire 4D space. The input kernel K(q, q') in BSE is taken as 
one-gluon-exchange like as regards color [(A^/2) • (A^ /2)] and spin (7^7^ ) dependence. 
The scalar function V(q — q') is a sum of one-gluon exchange Vqge and a confining term 



V con f. Thus we can write the interaction kernel as 



12|: 



K( q ,q>) = Qa«) ■ Qa<») Vf>V}*>V{Q-<n 



V(q-q') 



{q-q'Y + i W * /dr 



;i + 4a m 1 m 2 M 2 r 2 )- 1/2 - -§ 

UJn 



J(q—q')-r. 



12tt 



33 - 2f V A 2 



In 



Ml 



2\ "I 



M> = Max(M, m 1 + m 2 ). 



The Ansatz employed for the spring constant ui 2 ^ in Eq. (fT2l is [2j, [12J 



(12) 



uj 2 ^ = Arhi^M^lotsiM. 



(13) 



where rh\, m 2 are the Wightman-Garding definitions of masses of constituent quarks defined 
earlier. Here the proportionality of u 2 ^ on ots{Q 2 ) is needed to provide a more direct QCD 
motivation to confinement. This assumption further facilitates a flavour variation in u 2 ^. 
And uj 2 in Eq. ( fl~2i) and Eq. ( fl3l) is postulated as a universal spring constant which is 
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common to all flavours. Here in the expression for V(q — q'), as far as the integrand of the 
confining term V con f is concerned, the constant term Cq/uJq is designed to take account of 
the correct zero point energies, while a® term (a <C 1) simulates an effect of an almost 
linear confinement for heavy quark sectors (large m 1; m 2 ), while retaining the harmonic 
form for light quark sectors (small m 1; m 2 ) 12j as is believed to be true for QCD. Hence the 
term r 2 (l + 4aomim 2 M 2 ,r 2 ) _1//2 in the above expression is responsible for effecting a smooth 
transition from harmonic (qq) to linear (QQ) confinement. The basic input parameters in 
the kernel are just four i.e. ao = 0.028, Cq = 0.29, u = 0.158 GeV and QCD length scale 
A = 0.20 GeV and quark masses, m u ^ = 0.265 GeV, m s = 0.415 GeV, m c = 1.530 GeV and 



m b = 4.900 GeV which have been earlier fit to the mass spectrum of qq mesonsll2| obtained 
by solving the 3D BSE under Null-Plane Ansatz (NPA). However due to the fact that the 
3D BSE under CIA has a structure which is formally equivalent to the 3D BSE under NPA, 
near the surface P.q = 0, the qq mass spectral results in CIA formalism are exactly the 
same as the corresponding results under NPA formalism^, [l^j]. The details of BS model 
under CIA in respect of spectroscopy are thus directly taken over from NPA formalism (see 



Now comes to the problem of the 3D BS wave function. The ground state wave 



function <j)(q) satisfies the 3D BSE on the surface P ■ q — 0, which is appropriate for making 



contact with 0(3)-like mass spectrum (see 12j). Its fuller structure is reducible to that of 



a 3D harmonic oscillator with coefficients dependent on the hadron mass M and the total 
quantum number N. The ground state wave function <p(q) deducible from this equation thus 



has a gaussian structure [2|, ll2j and is expressible as: 

0(g) ~ e-?! 202 . (14) 

In the structure of 0(g) in (fT4"j) . the parameter f3 is the inverse range parameter which 
incorporates the content of BS dynamics and is dependent on the input kernel K(q, q'). The 

structure of the parameter j3 in <fi(q) is taken as 

o, {i n 

(3 2 = (2m 1 m 2 M^/ 7 2 )^ ; f = i _ (i 5 ) 

IVl >Co>q 

We now give the calculation of leptonic decays constants of pseudoscalar mesons em- 
ploying both LO and NLO Dirac covariants according to our power counting scheme in the 
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framework discussed in next section. 

III. CALCULATIONS AND RESULTS FOR f P 

A. Leptonic decays of pseudoscalar mesons to NLO 

Decay constants f P can be evaluated through the loop diagram which gives the coupling 
of the two-quark loop to the axial vector current and can be evaluated as: 

fpP^(0\Qij^ 5 Q\P{P)), (16) 

which can in turn be expressed as a loop integral, 

f P P, = J d 4 q Tr[^ P (P, <z)i 7M75 ]. (17) 

Bethe-Salpeter wave function ty(P,q) for a P-meson is expressed as, 

*(P,q) = S F (j> 1 )r(q)S F (-p 2 ), (18) 

which is expressed as the quark and anti-quark propagators flanking the Hadron-quark vertex 
T(q) function which is in turn expressed by Eq. (I91TT0]) . 

Using \&(P, q) from Eq. (ITS]) , and incorporating Hqq vertex function T(q) from Eq. (|9]IT0]) 
in Eq. (ITT]) , evaluating trace over the gamma matrices and multiplying both sides of Eq. 
( IT7j) by P M /(— M 2 )), we can express the leptonic decay constant f P as, 

f P = f!P + fP + f?+f?, (19) 
where f P , f P , f P , f P , are the contributions to f P from the four Dirac covariants asso- 
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ciated with coefficients Bi (i — 0, 1, 2, 3), and are expressed as: 



f (o) 
J p 



= V3N P B Jd 3 qD{q)<f>(q) f 



Mda 
2niA 1 A 2 



in 



+2j^ + A(m 1 -m 2 )a 



Tfti 

-2m x + 2j± - 2m 2 



m\m2 mim2 



M 2 



M 2 



f (i) 
J p 



V3N P B 1 J d 3 qD(q)(f>(q) J 
4 

+ {ml - m 2 )a 4Ma 

v 2 M 



Mda 
2mA 1 A 2 



M 3 M M 3 



4 "2 

m\ _ ^q z 



M 3 



M 



f (2) 



= V3N P B 2 J(PqD(q)c(>(q) / 



, . , mf mimi m 2 m 2 



2mA x A 2 

2 
2 



M 3 

4 



/ 2 2\ -2 



m 2 \ i /i *? 2 2,^/2 2\ 2 

-— nr <j + 4— — o" + — (m 2 - m> 

M 3 / M M 2 17 



/(f) = V3N P B 3 fd 3 qD(q)<f>(q) j 



Mda 
2-KiA 1 A 2 



mi + m2 .2 



(20) 



In deriving the above expressions, we had made use of the scalar products of various momenta 
expressed in terms of integration variables q and a as, 



Pi -P2 
Pi-P 

P2-P 

P-q 

2 
Pi 

Pi 

pi ■ q 

P2 ■ q 



-M 2 {m x + a)(m 2 -a)- q 3 , 
-M 2 (m 1 + a), 
-M 2 (m 2 -a), 
-M 2 a, 

-M 2 (mi + a) 2 + q 2 , 
-M 2 (m 2 -a) 2 + q 2 , 
1 {2g 2 - a\m\ - m\ + M 2 (l + 2a)]}, 
{-2q 2 + a[mj - m\ + M 2 (-l + 2a)]}. 



(21) 



We see that on the right hand side of the expression for f P , each of the expressions mul- 
tiplying the constant parameters B and B\ consist of two parts, of which only the second 
part explicitly involves the off-shell parameter a. It is can be seen that the off-shell con- 
tribution which vanishes for m\ = m 2 in case of using only the leading covariant 75, would 
no longer vanish for m\ = m 2 in the above calculation for f P (when other covariants are 
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incorporated in Hqq vertex function besides the leading covariant 75) due to the terms like 
4M and 4q 2 /M multiplying a 2 in fjp and fp^ respectively. This possibly implies that when 
other covariants besides 75 are incorporated into the vertex function, the off-shell part of 
fp does not arise from unequal mass kinematics alone (which is in complete contrast to the 
earlier CIA calculation of fp employing only 75). This may be a pointer to the fact that 
Dirac covariants other than 75 might also be important for the study of processes involving 
large q 2 (off-shell). Carrying out integration over da by method of contour integration by 
noting the pole positions in the complex a-plane: 



a -I- 0J\ ^ 9 9 -s9 

Ai = =>• a 1 = ±— -m^je, lu 1 = m 1 + q , 



x ^2 9 9 9 

=>• crj = Tjj +m 2 ±ie, u 2 = m 2 + q , 



(22) 
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we can again express f P as f P = f p 0) + /« + f p 2) + f p 3 \ where now 



fP = VZNpB J d 3 qD(q)<P(q) 



m? m 2 m 2 mvmi 

-2mi + 2— ^ - 2m 2 - 2^—^ - 2^—^ 

M 2 M 2 M 2 



+2 W 2 )W) +4{mi ' m2)Rl 



/£> = ^ NpBl J d 3 qD(q)<f)(q) 



M 3 M 



m\m 2 m?m| m 



+ 2^^ 1 



M 3 M 3 y £>($) 



f (2) 



V^JVp-Bi / d 3 qD(q)j>(q) 



4^2 2\ -2 1 



(23) 



/ ml m^m 2 mini- m~ , 



M 3 



M 



M 3 M 3 



+4—R 1 + — (m 2 . - m 2 )R 2 
M M V 2 



/<?> = v^iVp^g / d 3 qD(q)<f>(q) 



J_( mi + ma ) ff _j 



and -D((?) is given in Eq. ([6]), and the results of a-integration in the complex cx-plane, on 
whether the contour is closed from above or below the real <r-axis is: 



Ri 



Mda 
2mSAiA 2 ' 



-a 



M 2 (-c^i + u 2 ) + (mj - mf)(ux + u 2 ) 
\M 2 u x uj 2 \M 2 - (lu 1 + uj 2 ) 2 ] 



Ri 



Mda 
27riAiA 2 



a 



— M 4 - m{ 2 5m 2 + 4M 2 uxU 2 ){uJi + w 2 ) + 2M 2 mi 2 5m(a; 2 - wi) 
8MV^ 2 [M 2 - (cji + cj 2 ) 2 ] 



(24) 



To calculate BS normalizer N P for a pseudoscalar meson in the expression for f P in Eq. 
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(|23|) . we use the current conservation condition [2], 

2iP, = (27r) 4 | A Tr *(P,q) f^Spfa)) V(Pq)S F 1 (-p 2 ) 



+ (1<^2). (25) 



Putting BS wave function \I/(P, q) from Eq. (|T8|) in the above equation, carrying out deriva- 
tives of inverse of propagators of constituent quarks with respect to total momentum of 
hadron P^, evaluating trace over the gamma matrices, following usual steps and multiply- 
ing both sides of equation by P M /(— M 2 ) to extract out the normalizer Np from the above 
expression, we then express the above expression in terms of integration variables q and 
a. Noting that the four dimensional volume element d A q = d 3 (jMda, we then perform pole 
integration over da in complex cx-plane, making use of the pole positions in Eq. fT22|) . The 
calculation of normalizer is extremely complex due to unequal mass kinematics. We thus 
give here a general expression for the normalizer integral of the form, 

Np 1 = -(2tt) 2 * J d 3 qD 2 (g)0 2 (g)[^i(S, q)Ii + g 2 (B, q)I 2 + g 3 (B, q)I 3 +g 4 (B, q)h), (26) 

where B = (B , Bi, B 2 , B 3 ) and gi,...g^ are extremely complicated functions of B and q 
and are extremely lengthy expressions, and hence we do not present their actual forms here, 
whereas I\, ...I4 are analytic results of pole integration over the off-shell variable a in the 
complex cr-plane and are expressed as: 



h 



+00 



Mda 



2-ni 



2ui\ — M 2 ui 2 + 5ufuj 2 + 4co>iCo>2 + w| 
4^ 2 (M 2 -(^+u; 2 ) 2 ) 2 



Mda 
AfAl 



-a 



2ni- 



-M 4 uo 2 + (mj - ml)(uJi + uj 2 ) 2 (2uj 1 + u 2 )M 2 [6ujf + 9uo\uj 2 + 4wio;| + uj 2 (-m{ + mj + uo%)} 



8M 2 ufuj 2 [M 2 - (ut + u 2 ) 



212 
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+oo 

h = J 

— oo 

+M 4 [2co>i - 2m\u 2 + 2m\uj 2 + u\uj 2 + Auj x uj 2 2 + lu^} 

-M 2 [m\u 2 + m 2 u 2 + Auj\bJ 2 {u) x + uj 2 ) 2 + 2m 2 2 {-2u\ + io\w 2 + Au x u 2 2 + wf) 
—2m x {—2<jj\ + m 2 ci;2 + u x <jj 2 + Aujiuj 2 + u 2 )\} 



Mdo „ 



2?ri 



A 2 A 2 

(M 2 - m? + m 2 , + 2Mu 2 y 



8M 6 u 2 (M 2 - cu 2 + 2Mcu 2 + u 2 ) 2 

[M 2 + m\-m\- 2Mu l ) 2 [M i + M 2 (m\ -m\-J\- u 2 ) + (m 2 - m 2 )(3u; 2 - c^ 2 

16M 6 uf(M 2 - 2Mu 1 + cj 2 - u 2 2 ) 2 

(M 2 + m\-m 2 2 - 2Mu 1 ) 2 [-AMuji(m\ - m 2 2 + o;|)] 
16M 6 c^(M 2 - 2Mcji + cu 2 - cj 2 ) 2 



(27) 

After this, numerical integration over the 3-D variable d 3 q in Eq. (1261) is performed to 
evaluate Np. 

We have thus evaluated the expressions for fp and Np in framework of BSE under CIA, 
with Dirac structures of eq. ({TO]) , introduced in the Hqq vertex function besides 75 according 
to our power counting rule. We see that so far the results are independent of any model 
for <p(q). However, for calculating the numerical values of these decay constants one needs 
to know the constant coefficients Bo, B\, B 2 , B3 which are associated with the above Dirac 
structures. Because of the normalization condition, we take Bq = 1, and then there are 3 
parameters Bi/B , B 2 /B , B 3 /B Q , which will still be denoted as Bi, B 2 , B 3 for simplicity. To 
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see the contribution of various Dirac covariants on the calculation of meson decay constants, 
we first discuss the numerical procedure adopted to fit these coefficients 



B. Numerical Calculation 

Eq. (|23|) which expresses decay constants fp of pseudo-scalar mesons in terms of the 
parameters B , Bi, B 2 , B 3 is a highly non linear function of the B^s. This obviously implies 
that numerical methods must be applied to solve the problem. 

We used a simple Mathematica procedure for calculating the numerical integrals and 
searching for accurate values of the Bi (i = 0, ...3). We defined the following auxiliary 
function W(B) which is positive definite as, 

W(B) = J2lfp(B) ~ fp(exp.)} 2 , (28) 
p 

where B = (Bq, Bi, B 2 , B 3 ), and summation in the above equation runs over five pseu- 
doscalar mesons ir, K, D, Ds and B mesons studied in this work, and fp(exp.) are the 
central values of experimental data on decay constants 3 (indicated in Table II). 

From the numerical point of view the problem reduces to finding values of B^s such that 
W(B) has a minimum. We used Mathematica package which has some useful functions for 
minimizing. Those functions start from a point and search for a minimum near to that 
initial point. We constrained all the Bi's to lie within the interval [0,1]. We generated in a 
random way values of the Bi in this interval. Starting from those values, the Mathematica 
minimization function finds a minimum. Then it is checked if this minimum is "sufficiently 
near to zero" . This check is done by evaluating the percent average of the absolute values 
of the differences between the predicted fp values from the experimental value fp(exp.). 
Using this method we found that the values of coefficients B ,...,B 3 (with average error 
with respect to the experimental data less than 3.5%) respectively are: B$ = 1, B\/Bq = 
0.3727, B 2 /B = 0.2234, B 3 /B = 0.0821 to give the decay constant values, f v = 0.130 
GeV, f K = 0.164 GeV, f D = 0.194 GeV, f Ds = 0.296 GeV. and f B = 0.228 GeV which 
are within the error bars of experimental data Q] depicted in Table II for these five 
pseudoscalar mesons. These values of fp along with the contributions from various covariants 
and comparison with various models and experimental results are listed in Tables I and II. 
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There is one important point which needs to be clarified: The experimental data have 
different error bar, e.g., the data of n has a very high precision to the order of 0.1%, while 
for the case of B, the relative error is more than 16%. So in the fitting, we should take into 
account the difference, e.g., assign different weight for these data. However, we only give 
our formulation at NLO. From the above discussions, it is straightforward to rcognize, the 
smaller the meson mass, the larger the contributions of higher orders. For the case of pion, 
we even can perspect that higher order contributions (coming from higher order terms of 
Taylor series of B[s as powers of could be also very important. So, it is not reasonable 
to expect the NLO formulae can fit the data of pion to the precision of order of 0.1%. This 
is the reason why we fit the central value of data equally, as described above. 

IV. RADIATIVE DECAY CONSTANTS OF NEUTRAL P MESONS 

In this section we calculate the radiative decays of a neutral pseudoscalar meson such as 
7T° or i] c proceeding through the process P — ► 77 which proceed through the famous quark- 
triangle diagrams in the above framework using both the leading order and the next-to- 
leading order covariants in the Hadron-quark vertex function, taking the values of parameters 
B = 1, BJBq = 0.3727, B 2 /B = 0.2234, B 3 /B Q = 0.0821 fixed above in the calculation of 
fp values of n, K, D, D s and B mesons. The invariant amplitude for the decay of a neutral 
P-meson into two photons can be expressed as summation over the two triangle diagrams 
corresponding to the Direct and Exchange processes as: 



A(P^2 7 ) = 4= / d 4 qTr[^(P,q)i 1 .e 1 S F (q-Q)i 1 .e 2 }+-^= / d 4 qTr[V(P, q)^ .e 2 S piq+Q)^ .e± 



explicitly in Eq.(20) and Eq.(12)- (13), S F (q ± Q) are the propagators of the third quark 
in the Direct and Exchange diagrams respectively, where Q = k\ — k 2 is the the difference 
in momenta of the two emitted photons with momenta fciand k 2 respectively , while e\ t2 
are the polarization vectors of the two emitted photons in the above diagrams which differ 
from each other in the interchange 1 4=> 2. Evaluating traces over the gamma- matrices, 
combining various terms and then performing pole-integrations in the complex a— plane, we 





(29) 



where *&(P, q) is the BS wave function of a neutral P-meson given 
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can express amplitude for the above process as: 



A(P -> 2 7 ) = [F P }e^ pa P^ 2u Q p e la , (30) 

where P = pi + p 2 is the total hadron momentum, where pi,2 are the momenta of the 
quarks constituting the hadron, and the radiative decay constant, Fp is given as (the B 2 
term vanishes because of wrong charge parity), 



e 2 N F 



— 16m 4 4 8m 

B [8mS 1 ] + B 1 [- jr S 1 + -S 2 + -S 3 ] + ^[^(^ + S 3 + S 4 

(31) 

where 5*1,2,3,4,5 are the analytical results of integrals 



over the off-shell parameter a: 



s = f Mda = 12 (32) 

1 27riA 1 A 2 A 3 M 4 cj-20M 2 cu 3 + 64w 5 ' 1 ; 

(33) 



-J. 
-L 

-I 



00 
+00 



Mda 4 



^ 27riA 2 A 3 -M 2 u + 16cu 3 ' 
+°° Mrfa 4 



2 7 riA 1 A 3 -M 2 cj + 16c; 3 ' 



00 

+ °° Mda _ -1 

2vriA 1 A 3 (7 ~ -M 2 cu+ 16cj 3 ( ^ 

^ 5 " L 2«Aa' " -M 2 c + 16^3 {ib) 



evaluated by the method of contour integrations by noting the various pole positions in 
the complex a— plane : 
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& l = 0^a± = ±^- l -Tie; (36) 

A 2 = => at = + - =F ie; 

A 3 = =► af = T ie; uj 2 = m 2 + f 

corresponding to inverse propagators of the three quarks (of which A 12 correspond to 
the two constituent quarks in the meson) in the quark-triangle diagrams, expressed in terms 
of the off-shell parameter eras: 



Ax =LU 2 -M 2 (^ + a); (37) 
A 2 = u 2 -M 2 (±-a); 
A 3 = lu 2 - M 2 a 2 



From Eq.(33), it can be noticed that the contribution to radiative decay constant Fp 
from one of the next-to-leading order covariants associated with the parameter incompletely 
vanishes after trace evaluation. Numerical evaluation of Fp for tt° and rj c using the same 
set of parameters, E>i / B fixed from the calculation of leptonic decay constant fp values of 
7r, K, D, Ds and B mesons above gives F n = .031GeV -1 , F Vc = .006GeV _1 . These are very 
close to the experimental numbers F n (Exp.) = .025GeV _1 , and F Vc (Exp.) = .OOTAGeV' 1 
which are arrived at through the expression, T(P — > 27) = connecting the decay width 
T with radiative decay constants, Fp, using the central values of experimental data on decay 



widths for it and r) c mesons as r(7r° — > 27) = 8.5eV and T(?7 C — > 27) = 7AKeV 17|, Il8 ] 
respectively. 
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V. DISCUSSION 



In this paper we have calculated the decay constants fp of pseudoscalar mesons tt, K, 
D, Ds and B and radiative decay constants Fp for neutral pseudoscalar mesons tt° and 
r] c proceeding through the process P — > 27 in BSE under CIA. The Hadron-quark vertex 
function incorporates various Dirac covariants order-by-order in powers of inverse of meson 
mass within its structure in accordance with a power counting rule from their complete set. 
This power counting rule suggests that the maximum contribution to any meson observable 
should come from Dirac structures associated with Leading order terms alone, followed 
by Dirac structures associated with Next-to-Leading Order terms in the vertex function. 



fp values much closer to 



Incorporation of all these covariants is found to bring calculated 

results of experimental data Q and some recent calculations 5, (j, [l5|, [l6| for tt, K, D, 
Ds and B mesons. The fp are within the error bars of experimental data for each one of 
these five mesons by fitting three parameters. The calculation of radiative decay constants 
of 7T° and 1] C is again close to the experimental data[l3, Q- 

The results for tt, K, D, Ds and B mesons with parameter set: B = 1, Bi/B = 0.3727, 
B 2 /B Q = 0.2234, B 3 /B = 0.0821 (giving f P values with average error with respect to 
experimental data less than 3.5%) are presented in Table I. In Fig. 1 we are plotting functions 
I P (q) (i = 0, ...3) vs q, where P P (q) is the integrand of fp in equations fl23l . The plots of 
variations of I p (q), ...Ip(q) with q for tt, K, D, Ds and B mesons, along with the results 
in Table I, show that the contribution to fp from NLO covariants is much smaller than the 
contribution from LO covariants for K, D, Ds and B mesons. Comparison with experimental 
data and other models is shown in Table II. It is seen from Table I that as far as the various 
contributions to decay constants fp are concerned, for K mesons, the LO terms contribute 
60%, while NLO terms 40%. However for heavy-light meson D, the LO contribution increases 
to 90%, while NLO contribution is 10%. For Ds meson, LO contribution is 91%, while NLO 
contribution is 9%. But for B meson, the LO contribution is 96%, while NLO contribution 
reduces to just 4%. This is in conformity with the power counting rule according to which 
the leading order covariants, 75 and 275 (7 • P)(l/M) (associated with coefficients B and 
Bi) should contribute maximum to decay constants followed by the next-to-leading order 
covariants, —^75(7 • q){l/M) and — 75 [(7 • P)(l ' <?) — (7 ' Q)(l ' P )]0-/M 2 ) (associated with 
coefficients B 2 and -B3) in the BS wave function, Eq. (j9l-f fT0l) . 
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However the situation is different for the lightest meson tc which enjoys a unique status 



due to the fact that mass of a pion, M is unusually small (<< Aqcd), and the large difference 
between the sum of two constituent quark masses and the pion mass shows that the quarks 
are far off shell and the internal momentum q should be the same order as the pion mass 
and the approximation q « P ~ M breaks down for pion. Hence the contribution of NLO 
covariants in pion case is even larger than the contribution of LO covariants. Thus, the 
NLO covariants in pion should play a more dominant role in contrast to heavier mesons K, 
D, D$ and B. However the sum of LO and NLO contributions adds up to the experimental 
value for pion fp (=0.130 GeV). Further investigations on higher order terms can show even 
more details of the pion structure. 

To check the validity of our calculation, we then do numerical evaluation of radiative 
decay constants Fp for tt° and r\ c using the same set of parameters, Bi/Bofixed above from 
the calculation of leptonic decay constant fp values of tt, K, D, D$ and B mesons. This 
gives F n = .031GeV~ l , F Vc = .006GeV _1 . These are very close to the experimental numbers 
F v (Exp.) = .025GeV~ 1 , and F Vc (Exp.) = .0074Gel / ~ 1 which are arrived at through the 
expression, T(P — > 27) = -fj— connecting the decay width T with radiative decay constants, 
Fp, using the central values of experimental data on decay widths for ir and rj c mesons as 
r(7r° -> 2 7 ) = 8.5eV and T(t] c -> 27) = 7.4 [l7, Q respectively. 

The numerical results for leptonic decay constants, fp and radiative decay constants, Fp 
obtained in our framework upto the next to leading order covariants demonstrates the valid- 
ity of our power counting rule, which also provides a practical means of incorporating various 
Dirac covariants in the BS wave function of a hadron. By this rule, we also get to understand 
the relative importance of various covariants to calculate various meson observables. This 
would in turn help in obtaining a better understanding of the hadron structure. Here would 
would like mention the robustness of our framework: On one hand, at lower order (s), with 
limited number of parameters, we can globally reproduce almost all the decay constants of 
certain kinds of meson. On the other hand, by introducing higher order corrections, we can 
accommodate enough parameters to fit the data as precise as possible, so than to get a good 
parameterization of the structure of certain special hadron for further investigations. 
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fp 


fp 


\f L P°\ 


\fr°\ 






-p _ fLO i f NLO 
JP — Jp + Jp 


71 


0.110 


-0.154 


0.000 


0.175 


0.044 


0.175 


25% 


75% 


0.130 


K 


0.202 


-0.104 


0.025 


0.039 


0.098 


0.064 


60% 


40% 


0.164 


D 


0.271 


-0.097 


0.010 


0.009 


0.174 


0.019 


90% 


10% 


0.194 


D s 


0.426 


-0.156 


0.013 


0.013 


0.270 


0.026 


91% 


9% 


0.296 


B 


0.345 


-0.125 


0.005 


0.003 


0.220 


0.008 


96% 


4% 


0.228 



TABLE I: Decay constant fp values (in GeV) for ir ,K ,D Ds and B mesons in BSE with the 
individual contributions fp, fp, fp, fp from various Dirac covariants along with the contributions 
from LO and NLO covariants and also their % contributions for parameter set: -Bo = L B\/Bq = 
0.3727, B2/B0 = 0.2234, B^/Bq = 0.0821 (with average error with respect to the experimental 
data less than 3.5%) 
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fx 


fo 


fn s 


f B 


BSE (3.5% average error) 












present paper 


0.130 


0.164 


0.194 


0.296 


0.228 


BSE [5] 








0.248 




SDE [6] 




0.164 








Lattice [15] 






0.208±0.004 


0.241±0.003 




QCD-SR [16] 






0.20±0.02 


0.23±0.02 




Exp. Results [13] 


0.1300±0.0001 


0.159±0.001 


0.22±0.02 


0.29±0.03 




Babar+Belle 












Collaboration [14] 










0.24±0.04 



TABLE II: Comparison of results of fp (in GeV) for tx K, D, Ds and B in BSE with the parameter 
set B = 0.7045, B x = 0.2626, B 2 = 0.1574, B 3 = 0.0579 (with average error 3.5%) with those of 
other models and experimental data. 
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